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Abstract. We report on the first observation of the 
Seyfert 1 galaxy NGC 4593 in the 0.1-200 keV band, 
performed with the BeppoSAX observatory. Its spectral 
components are for the first time simultaneously mea- 
sured: a power-law with photon spectral index T ~ 1.9; 
the Compton-reflection of the primary power-law; a mod- 
erately broad (p > 60 eV) K Q fluorescent line from neu- 
tral iron; and an absorption edge, whose threshold en- 
ergy is consistent with K-shell photoionization from Ovil. 
The amount of reflection and the iron line properties are 
consistent with both being produced in a plane-parallel, 
X-ray illuminated relativistic accretion disc surrounding 
the nuclear black hole, seen at an inclination of ~ 30°. 
Any cutoff of the intrinsic continuum is constrained to lay 
above 150 keV. The claim for a strongly variable soft ex- 
cess is dismissed by our data and by a reanalysis of archival 
ASCA and ROSAT data. 
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1. Introduction 

It is widely accepted that the huge energy output of Ac- 
tive Galactic Nuclei (AGN) is due to the release of grav- 
itational energy by matter falling onto a supermassive 
(~ 10 6 -10 8 M Q ) black hole. It has been known since the 
beginning of the X-ray astronomy that the spectrum of 
Seyfert 1 galaxies above a few keV and up to a few tens of 
keV can be described at the 0-th order by a simple power- 
law, with typical photon index T ~ 1.5 — 2 (Mushotzky 
1984). Later on, an emission line from neutral or mildly 
ionized iron (Holt et al. 1980; Perola et al. 1986; Pounds et 
al. 1990) and a flattening of the spectrum above ~ 10 keV 
(Pounds et al. 1990; Piro et al. 1990) were discovered to be 
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common spectral features in the class as well. Both these 
features have been interpreted as the effect of reprocess- 
ing of the primary radiation by optically thick matter sur- 
rounding the nucleus (George & Fabian 1991; Matt et al. 
1991). The detection of the double-horned and redshifted 
profile of the iron line in a long-look ASCA observation 
of the archetypical Seyfert 1 MCG-6-30-15 (Tanaka et al. 
1995; Iwasawa et al. 1996), together with the general evi- 
dence that the iron lines are on the average broad in this 
class (Nandra et al. 1997), gave the first direct evidence 
that the reprocessing matter is located close to the black 
hole, probably in a Keplerian accretion disk. 

In soft X-rays, an excess above the extrapolation of 
the high-energy power-law was measured in 50% of the 
Seyferts observed by EXOSAT (Turner & Pounds 1989) 
and about 90% of those observed by ROSAT (Turner et 
al. 1993; Walter & Fink 1993). The good correlation with 
Optical/UV (Walter & Fink 1993; Puchnarewicz et al. 
1996; Laor et al. 1997) has traditionally supported the 
idea that this excess represents the hard tail of thermal 
emission from the accretion disk (Czerny & Elvis 1987), 
although in at least a few cases reprocessing by ionized 
matter seems a more viable explanation (Piro et al. 1997; 
Guainazzi et al. 1998b). However, caution must be cm- 
ployed when evaluating the soft excess in band-limited 
detectors and/or from multi-instrumental fits of not si- 
multaneous observations. The extrapolation of 1-10 keV 
power-law fits can yield "faked" soft excesses, because the 
so determined spectral indices can be harder than the true 
one, due to the flattening contribution of the reflection 
component. Eventually, the discovery of absorption edges 
from highly ionized species of oxygen in almost 50% of 
the Seyfert Is observed by ASCA so far (Reynolds 1997; 
George et al. 1998) has suggested that in some cases the 
apparent soft excess could be simply due to a mis-fit of a 
complex and ionized absorber. 

Given the spectral complexity outlined above a de- 
tailed and self-consistent description of the X-ray spec- 
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tra of Seyfert Is requires broadband spectral observations. 
Moreover, the various spectral components are expected to 
be produced in different physical regions (see Mushotzky 
et al. 1993 for a review of the pre-ASCA interpretative 
scenario, which is still largely valid) and are therefore ex- 
pected to exhibit different variability patterns and/or a 
delayed response to the changes of the nuclear radiation 
flux. The simultaneous measure of all the spectral compo- 
nents is a crucial requirement for any X-ray observation 
of Seyfert galaxies (see the discussion in Cappi et al. 1996 
about the dependence of the iron line measurements on the 
underlying continuum determination). The Italian-Dutch 
satellite BeppoSAX (Boclla et al. 1997a) carries a scien- 
tific payload which covers the unprecedented wide energy 
band between 0.1 and 200 keV, with imaging capabilities 
and good energy resolution in the 0.1-10 keV band. A pro- 
gram of spectral survey of a sizeable sample of Seyfert 1 
galaxies is ongoing, and in this paper we report the results 
of the observation of NGC 4593. 

NGC 4593 (a 20 oo = 12 h 39 m 39 s .4, S 2 qoq = -5°20'39") 
is a nearby (z = 0.009) barred spiral galaxy of Hubble type 
SBb, which hosts a Seyfert 1 nucleus. The Spectral En- 
ergy Distribution (SED) is characterized by a very weak 
or missing "blue bump" (Santos-Lleo et al. 1994). In X- 
rays, the source is variable, on timescales of weeks- months, 
by a factor of ~ 3 in the intermediate (i.e.: 2-10 keV) 
and ~ 4-5 in the soft X-rays (Ghosh & Soundararajape- 
rumal 1993). In NGC 4593 a soft excess above the ex- 
trapolation of the intermediate X-ray power-law was ob- 
served both by EXOSAT (Ghosh & Soundararajaperu- 
mal 1993; Santos-Lleo at al. 1995) and ROSAT (Walter 
& Fink 1993). It was claimed to be variable, ranging from 
to ~ 270% of the intermediate X-ray extrapolated flux. 
However, this picture is further complicated by the ASCA 
discovery of a warm absorber (Reynolds 1997), with op- 
tical depths of the Ovil and Ovin photoionization edges 
equal to ~ 0.3 and 0.1, respectively. There was no signifi- 
cant evidence of broadening of a weak (Equivalent Width 
EW ~ 90 eV) fluorescent line from neutral iron (centroid 
energy E c ~ 6.35 keV) in the ASCA data (Nandra et al. 
1997). In the hard X-ray regime, BATSE measured a 20- 
100 keV flux of (9.8±1.6) x lO" 11 erg cm" 2 s" 1 (Malizia et 
al. 1997), whereas the 50-150 keV OSSE 2-cr upper limit 
was 4.4 x 10~ n erg cm~ 2 s _1 (Johnson et al. 1993). 

This paper is organized as follows. In Sect. 2 we de- 
scribe the data reduction and preparation. Sect. 3 and 4 
deal with the timing and spectral analysis of BeppoSAX 
data, respectively. In Sect. 5, we compare our findings with 
a reanalysis of archival ASCA data of the same object. The 
results are discussed in Sect. 6. 

2. Data Reduction 

The Italian-Dutch satellite BeppoSAX carries four co- 
aligned Narrow Field Instruments. Two of them are gas 
scintillation proportional counters with imaging capabili- 



ties: the Low Energy Concentrator Spectrometer (LECS, 
0.1-10 keV, Parmar et al. 1997) and the Medium Energy 
Concentrator Spectrometer (MECS, 1.8-10.5 keV, Boclla 
et al. 1997b). The other two instruments arc direct- view 
detectors, seen through a rocking collimator to achieve a 
continuous monitoring of the background: the High Pres- 
sure Gas Scintillator Proportional Counter (HPGSPC, 4- 
120 keV, Manzo et al. 1997) and the Phoswitch Detec- 
tor System (PDS, 13-200 keV, Frontera et al. 1997). The 
HPGSPC is tuned for spectroscopy of very bright sources 
with good energy resolution, while the PDS possesses an 
unprecedented sensitivity in its energy bandpass. Only 
LECS, MECS and PDS data will be considered in this pa- 
per. The HPGSPC data points could not in fact provide 
significant constraints on the spectral shape of NGC 4593. 

NGC4593 was observed by BeppoSAX from 1997 De- 
cember 31 05:30:18 UT to 1998 January 2 07:10:17 UT. 
Data were telemetred in direct modes for all instruments, 
which provide full information about the arrival time, en- 
ergy, burst length/rise time (RT) and, when available, po- 
sition for each photon. Standard reduction procedures and 
screening criteria have been adopted to produce linearized 
and equalized event files. In particular, time intervals have 
been excluded from the scientific product accumulation 
when: the angle between the pointing direction and the 
Earth's limb was < 5°; the momentum associated to Ge- 
omagnetic Cutoff Rigidity was > 6 GcV/c; the satellite 
passed through the South Atlantic Geomagnetic Anomaly 
(SAGA). The PDS data have been further screened by 
eliminating 5 minutes after any SAGA passage to avoid 
gain instabilities due to the recovery to the nominal volt- 
age value after instrumental switch-on. The RT selection 
has been performed using crystal temperature dependent 
thresholds (instead of the fixed thresholds in the standard 
processing). Total exposure times were 31426 s, 84195 s, 
and 38467 s for the LECS, MECS and PDS, respectively. 

Spectra of the imaging instruments have been ex- 
tracted from circular regions of radius 8' and 6' around 
the apparent centroid of the source for the LECS and 
MECS, respectively. The background subtraction has been 
performed using spectra extracted from blank sky event 
files in the same region as the source. The background 
contributes 5% and 3% to the LECS and MECS count 
rates in the relevant energy band passes. PDS background- 
subtracted products (spectra and light curves) have been 
produced by plain subtraction of the "off-" from the "on- 
source" products. Total net count rates were 0.310 ± 
0.003 s-\ 0.350 ± 0.002 s" 1 and 0.80 ± 0.04 s" 1 in the 
LECS, MECS and PDS, respectively. 

In this paper: errors are quoted at 90% level of 
confidence for two interesting parameters (A% 2 = 4.61, 
Lampton et al. 1976); energies arc in the source rest 
frame; the cosmological parameters q = 0.5 and H = 
50 km s _1 Mpc" 1 are assumed, unless otherwise speci- 
fied. 
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Fig. 1. Light curves in the 0.7-3 keV (LECS, upper panel) 
and 3-10 keV (MECS, central panel) energy bands and 
the corresponding HR (lower panel). The binning time is 
2048 s. The light curves are not background-subtracted 
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Fig. 2. Light curves in the 3-10 keV (MECS. upper 
panel), 13-200 keV (PDS, central panel) energy bands and 
the corresponding HR (lower panel). The binning time 
is 25000 s. The MECS light curve is not background- 
subtracted 



3. Timing analysis 

NGC 4593 is well known to display rapid fluctuations in X- 
rays, with two-folding time ~ 1.1 h (Barr et al. 1987). The 
count rates in the BeppoSAX instruments span indeed a 
dynamical range of about a factor of 2 on timescales as 
low as 3 x 10 4 s. We have therefore extensively searched 
for spectral variability associated with the flux changes, 
extracting light curves and Hardness Ratios (HR) in dif- 
ferent energy bands and with various binning times. In 
Fig. [l] and we show the light curves in the 0.7-3 keV, 
3-10 keV and 13-200 keV with the corresponding contigu- 
ous HR. The bands have been chosen in order to sample 
energy intervals where different spectral components dom- 
inate, see Sect. 4. In particular, the softer band is affected 
by the warm absorber, which is basically transparent for 



Fig. 3. 0.1-3 vs. 3-10 keV HR versus the 3-10 keV count 
rate. Each data point corresponds to a 2048 s light curve 
bin 



photons of energy higher than 3 keV. No significant vari- 
ability of the HR is detected. If a fit with a constant is 
performed on the 3-10 keV vs. 0.7-3 keV HR, the x 2 is 
29.3 for 38 degrees of freedom (dof). The same quantity 
for the 13-200 keV vs. 3-10 keV HR is 4.8/6 dof. A tiny 
perturbation of the intermediate X-ray flux in a 2 x 10 4 s 
interval around 9 x 10 4 s after the beginning of the obser- 
vation is not accompanied by a similar phenomenology in 
the soft X-rays and yields therefore a decrease of the HR 
by a factor ~ 15%. However, no variation in the spectral 
parameters above the statistical uncertainties has been de- 
tected in time-resolved spectroscopy (assuming the base- 
line model of Sect. 4.1). Similar results are obtained if the 
bands 0.1-3 and 3-10 keV are considered (see Fig. ||). We 
will therefore focus in the following on the time-averaged 
spectral behavior. 



4. Spectral analysis 

The spectra of the imaging instruments have been re- 
binned in order to sample the intrinsic energy resolution 
of the detectors with 3 (LECS) or 4 (MECS) energy chan- 
nels. Each channel has at least 30 counts, which ensures 
the applicability of x 2 test. The PDS spectrum has been 
quasi-logarithmically rcbinned, in order to have 16 energy 
channels in the 14-200 keV band. The spectra of the three 
detectors have been fitted simultaneously. Numerical rel- 
ative normalization factors among the BeppoSAX instru- 
ments have been added to all the following spectral fits. 
The reasons is two-fold: a) the BeppoSAX instrument re- 
sponse matrices employed in this paper (September 1997 
release) exhibit slight mismatches in the absolute flux cal- 
ibration; b) the sampling of the instruments is not strictly 
simultaneous, due to the need for operating the LECS only 
during satellite nights or the different data selection crite- 
ria between imaging instruments and the PDS. This can 
affect the flux measured in variable sources as NGC 4593. 
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Fig. 4. Spectra and best-fit model (upper panel) and 
residuals in units of standard deviations (lower panel) 
when a simple power-law model with photoelectric ab- 
sorption is applied 

The LECS to MECS factor has been left free in the fit- 
ting procedure, and turns out to be comprised in the range 
0.74-0.76, which is consistent with typical values observed 
so far (0.7-1.1: Grandi et al. 1997; Haardt et al. 1998; 
Cusumano et al. 1998). The PDS to MECS factor has 
been instead held fixed to 0.8, as the available statistics 
was not good enough to provide independent constraints 
on it. This value corresponds to the multiplication of the 
best fiducial value estimated by Cusumano et al. (1998) 
by 0.82, to account for the effect of the PDS RT selection 
algorithm employed (see Sect. 1). The systematic uncer- 
tainty on this parameter can be estimated ±10%. Spectral 
fits have been performed in the 0.1-4 keV (LECS), 1.8- 
10.5 (MECS), 14-200 keV (PDS) energy bands. 

4-1- Continuum shape 

In Fig. ^ the result is shown, when a simple power-law 
model with photoelectric absorption is applied. The qual- 
ity of the fit is rather poor (\ 2 = 262.6/122 dof). The 
main deviations are due to: (a) an absorption feature 
starting at E ~ 0.7 keV; (b) a prominent emission line 
with centroid energy E ~ 6.5 keV; (c) a "bump" in the 
PDS band, peaking at energy E ~ 25 keV. The residu- 
als between 0.3 and 0.6 keV (immediately red wards the 
absorption feature) are systematically positive. In princi- 
ple, this may be due to the emergence of a soft excess 
in this energy range. However, the recovery of the residu- 
als to values consistent with zero below 0.3 keV leads to 
the suggestion the the above feature is simply the typical 
wavy residual produced by a mis-fit absorption edge (see 
e.g. Nandra & Pounds 1992). The BeppoSAX broadband 
allows to study simultaneously the whole spectral com- 
plexities expected on the basis of previous band-limited 
measures of this object and of the Seyfert Is as a class (cf. 
Sect. 1). We have then defined a "baseline" model, where a 



Compton reflection component from a neutral slab (model 
pexrav in Xspec, Magdziarz & Zdziarski 1995), an emis- 
sion line and a photoionization absorption edge are su- 
perposed to the photoelectric absorbed power-law. The 
possibility that the reflecting matter is substantially ion- 
ized is not required by the data. The model depends on 
the heavy element abundance (which has been held fixed 
to the solar one), and on the angle between the normal 
to the slab and the line of sight ("inclination angle" 6). 
has been held fixed to 30° hereinafter (the best-fit value 
arising from a fit of the iron line profile with a relativistic 
model, see Sect. 4.2). The only free parameter in addition 
to the power-law model of the continuum is the relative 
normalization R between the reflected and the primary 
components (equal to 1 for an isotropic source, illuminat- 
ing a plane-parallel infinite slab). The best-fit R and V 
are basically unaffected, if one assumes that the reflector 
is seen face-on (i.e.: 8 = 0°). The "baseline" model yields 
a very good x 2 = 99.0/116 dof. The Table | reports the 
best-fit parameters and results. 

It is worth noticing that no soft excess above the ex- 
trapolation of the high-energy power-law is required. A 
model, where one replaces the absorption edge with a con- 
tinuous break of the primary power-law provides a much 
worse fit (x 2 = 133.6/116 dof) and leaves significant resid- 
uals in the 0.2-1.3 keV energy range. Moreover, the best-fit 
spectrum is convex (i.e.: T so [ t < T), further demonstrat- 
ing that no soft excess is present in the data. The inferred 
absorbing column density is consistent with the Galactic 
contribution along the line of sight as measured by Elvis et 
al. (1989, Njj Gal = 1.97 x 10 20 cm" 2 ). The best-fit model 
and deconvolved spectra are shown in Fig. ^ The 0.1- 
2 keV, 2-10 keV and 20-100 keV fluxes are ~ 2.11 x KT 11 , 
3.74 x f0 -11 and 7.22 x f0~ n erg cm~ 2 s^ 1 , respectively. 
They correspond to rest frame luminosities of 5.80 x 10 42 , 
1.03 x 10 43 and 1.99 x 10 43 erg s _1 , respectively. 

The unprecedented BeppoSAX energy coverage allows 
the simultaneous determination of the primary radiation 
steepness and of the Compton reflection intensity with 
the best accuracy ever achieved. In Fig. || the contour 
plot for the photon index versus the relative normaliza- 
tion between the reprocessed and primary components R 
is shown. The latter parameter is 1 if the reflection occurs 
in an infinite, plane-parallel slab and the primary source 
emits isotropically. Higher values are formally possible, 
and may be due to the geometry of the reflecting mat- 
ter, covering more than 2tt (e.g. in a concave or warped 
accretion disk), to an intrinsic anisotropy of the primary 
source, or to a delayed response of the reflecting matter 
to changes of the primary flux. At 90% level of confidence 
for two interesting parameters, the photon index is con- 
strained between 1.81 and 1.91. R is consistent with a 
plane-parallel geometry of the reflecting matter (90% con- 
fidence interval between 0.6 and 1.6). Any cut-off of the 
intrinsic power-law is constrained to lay at energies higher 
then 150 keV (see Fig.@). It is worth noticing that the sta- 
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Table 1. BeppoSAX best-fit results, wa = photoelectric absorption from neutral matter; po = power-law; px - 
ton reflection; ga — Gaussian line; ed = photoionization absorption edge; bk = broken power-law 
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(keV) 


r or r soft 


X7 dof 
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1.55± ; 13 


1.70 ± 0.04 














262.6/122 


wa*px 


1.9 ± 0.2 


1.79 ±0.04 


o.9±8:l 












183.9/121 


wa* (px+ga) 


1.9 ±0.2 


1.80 ± 0.04 


o.8±8:l 


6.41±gif 


o.2±S!:l 


170±g8 






134.8/118 


wa*ed* (px+ga) a 


2.3 ±0.3 


1.87 ±0.05 


1.1 ± 0.4 


6.42±gil 


o.3±g:i 


190±|S 


0.77±°;°! 


0.43±H 6 7 


98.5/116 


wa*ed* (px+ga) 


2.3 ±0.3 


1.86 ± 0.05 


i-o±8:S 


6.42±°:i! 


0.3±8;| 


200±r 


0.78 ±0.06 


o.42±8:i6 


99.0/116 


wa*ed* (px+ga) c 


i-5±i:l 


i.96±8:8§ 


o.6±J:8 


6.5±8;1 


o.5±J:S 


280±fi8 


0.68 ± 0.03 

o.83±8:81 


o.30±8:os 
o.i3±8:8 5 7 


1396.4/1357 


wa* (bk+px+ga) 


1.4 ±0.3 


i.sT±8:oi 


i-2±S:l 


6.42±8i? 


o.4±8:§ 


240±53° 


2.2 ±0.5 


1.67 ±0.07 


133.6/116 



a 9 = 0; 

b 6 — 30°, "baseline model" in text. 
c fit on the ASCA data (see Sect. 5) 
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Fig. 5. Left: spectra and best-fit model (upper panel) and residuals in units of standard deviations (lower panel) 
when the "baseline" model is applied. Right: unfolded energy spectrum and best-fit model (solid lines). The direct 
and reflected continuum and the emission line are separately indicated with dotted lines. The location of the O vn 
photoionization absorption edge in the observers frame is labeled 



tistical relative uncertainties on T and R are rather small, 
despite the strong correlation between the two parame- 
ters (note the strongly inclined contour plot in Fig. Q). 
For comparison, the single parameter 68% statistical er- 
rors (Ar = 0.02 and AR/R = 0.17) are 4.5 and 3.5 times 
smaller than the ones obtained from the Ginga observa- 
tions of the same Seyfert 1 (Nandra & Pounds 1994). It 
should however be noticed that the residual systematic 
uncertainties on the relative PDS to MECS normaliza- 
tion factor affect the accuracy of both these parameters, 
with additional uncertainties of about 1% and 30%, re- 
spectively. 

4-2. On the K a iron line 

The centroid energy of the iron line is well consistent with 
Ka fluorescence from neutral iron. The line is broad if a 
simple broad Gaussian profile is used to describe it (see 
Fig. ||); the intrinsic width is comprised in the range GO- 
GOO eV at 90% level of confidence for two interesting pa- 



rameters (best- fit value ~ 300 eV). We have tried also 
alternative parameterization of the line. Adding a fur- 
ther "narrow" line (i.e.: intrinsic width a held fixed at 
zero) to the "baseline" model results in no improvement 
of the quality of the fit (A% 2 = 0). Given the intrinsic 
width of the line and the agreement between the derived 
best-fit line EW and the theoretical expectations if the 
line is produced in an X-ray illuminated relativistic ac- 
cretion disk (Matt et al. 1992), we have tried to use a 
self-consistent model of line emission from a relativistic 
accretion disk (model diskline in Xspec, Fabian et al. 
1989). If all its parameters are allowed to be free, most 
of them are totally unconstrained. No further constraint 
comes from the ASCA data as well (Nandra et al. 1997). 
If it is assumed that the inner radius Rj of the emitting 
region is 6 gravitational radii (Rq = GM/c 2 ); the emis- 
sivity law index is equal to -2.5 (Nandra et al. 1997); 
the inclination of the line and Compton reflection contin- 
uum emitting region is the same; and the line is neutral 
(i.e.: E c = 6.4 keV); then 6 = 32°± 2 3, log(R D ) > 2.1 and 
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Fig. 6. Contour plot of the relative normalization be- 
tween the primary and Compton reflected spectral compo- 
nents versus the intrinsic photon index, when the "base- 
line" model is employed. Iso-x 2 curves are at 68%, 90% 
and 99% confidence level for two interesting parameters 
(Ax 2 = 2.3, 4.6 and 9.2) 
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Fig. 7. Contour plot of the cutoff energy of the primary 
power-law component versus the photon index when the 
"baseline" model is employed 



EW = 240±go cV. The quality of the iron line modeling 
(x 2 = 98.9/1 f 4 dof) is comparably good as for the broad 
Gaussian. 

4-3. The warm absorber 

In 12 out of 24 Seyfert galaxies observed by ASCA, 
absorption edges from ionized species of oxygen have 
been detected (Reynolds 1997), which have been inter- 
preted as the imprinting of warm gas along the path 
from the nucleus and us. NGC 4593 is one of these ob- 
jects, and the BeppoSAX observation confirms this out- 
come, thanks to the detection of an absorption edge 
with threshold energy E th = 0.78 ± 0.06 keV. The edge 
energy is consistent with the K-photoionization thresh- 



Fig. 8. Contour plot of the intrinsic width versus the ccn- 
troid energy of the iron line, when the "baseline" model 
is employed 



old energy of Ovu. We have tried to give a qualitative 
characterization of the ionization and chemical structure 
of the absorbing matter, by tentatively including in the 
fit four absorption edges, with threshold energies held 
fixed to the values expected from Ovu (0.739 keV), Ovill 
(0.871 keV), Neix (1.196 keV) and Nex (1.362 keV). 
However, the available statistics is not good enough 
to give us significant constraints. Only the Ovu edge 
yields a significant detection (tovii = 0.32lg 'Jg), while 
only upper limits can be obtained for the other three 
edges (toviii < 0.22, r No i X < 0.05, r NcX < 0.03). Using 
Cloudy (Ferland 1996), we have constructed a grid of 
self-consistent models of the spectra transmitted through 
a ionized gas in thermal and ionization equilibrium, when 
the SED is the one observed in NGC 4593 (Santos-Lleo 
et al. 1995). They depend on the ionization parameter 
U (defined as the dimensionless ratio between the num- 
ber of Hydrogen ionizing photons and the electron den- 
sity of the gas), and the warm column density N\y- The 
fit is comparably good as the phcnomenological descrip- 
tion of the "baseline" model (x 2 = 104.6/116). The 
best-fit parameters of the warm absorbing matter are 
N w = (2.5±};|) x 10 21 cm" 2 and log(U) = -0.7 ± 0.3. 
For a source with a ox ~ 1.20, the best-fit U corresponds to 
a value of the more common £ = L/n 2 R ~ a few. The con- 
tinuum parameters are slightly affected (r = 1.92 ± 0.10; 
R = 1.3 ± 0.5), but remain consistent with the one of the 
"baseline" model within the statistical uncertainties. 



5. Comparison with ASCA and ROSAT results 

Nandra et al. (1997) detected with ASCA a weak 
(EW = 90±|q eV) and narrow iron line, whose best-fit pa- 
rameters are inconsistent with the ones corresponding to 
the "baseline" BeppoSAX best-fit model. We have there- 
fore re-analyzed the ASCA observation to check the ro- 
bustness of the iron line variability. ASCA (Tanaka et al. 
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Fig. 9. Spectra and best-fit model (upper panel) and 
residuals in units of standard deviations (lower panel) 
when the modified BeppoSAX "baseline" model is applied 
to the January 1994 NGC 4593 ASCA observation. Only 
SIS1 and GIS3 data are shown for clarity 

1994) payload include a pair of Charged-coupled devices 
(SISO and SIS1, 0.57-9 keV) and a pair of gas scintilla- 
tion proportional counters (GIS2 and GIS3, 0.7-10 keV). 
It observed NGC 4593 on January 9, 1994. Spectra have 
been extracted from the screened event files, which have 
been filtered according to standard criteria, using extrac- 
tion radii of 2.30', 3' and 6' for the SIS1, SISO and GIS, 
respectively. Net exposure times amount to 30 and 33 ks 
for the SIS and GIS, respectively. The spectra have been 
rcbinncd in order to have at least 20 counts per energy 
channel. Background spectra have been extracted from 
blank sky pointing event files, using the same extrac- 
tion region in detector coordinates as the source. The 
same BeppoSAX "baseline" model has been applied to 
the spectra of all detectors simultaneously, except for 
the fact that ASCA is capable to resolve the Ovu and 
Ovin absorption edges. The parameter R, which cannot 
be constrained by ASCA given the limited bandwidth 
of its instruments, has been forced to be comprised in 
the 90% confidence level range determined by BeppoSAX 
data, i.e.: 0.6-1.6. The time-averaged 2-10 keV flux in the 
ASCA observation is very close to that measured by Bep- 
poSAX (F ~ 3.70 x 10' 11 erg s' 1 cm' 2 ). No significant 
deviation exists between the ASCA and BeppoSAX best- 
fit parameters, with the exception of a slight steepening 
of the intrinsic power-law (AT = 0.10±§;j$, cf. Table @). 
In particular, the iron emission line best-fit parameters 
are consistent with the one obtained from the BeppoSAX 
data analysis within the statistical uncertainties and there 
is no evidence for a much weaker line in the ASCA data. 
It is worth noticing that no soft excess is required. The 
ASCA spectra and residuals against the best-fit are shown 
in Fig. |. 

Walter & Fink (1993) report the measure of a very 
steep (T = 2.5) intrinsic spectral index with ROSAT. 



However, their analysis does not take into account the pos- 
sible contribution of a ionized absorber. We therefore re- 
analyzed the data of a pointed NGC 4593 ROSAT/PSPC 
observation, to compare BeppoSAX and ROSAT find- 
ings. The ROSAT observation was performed July 14 
1992. The average spectrum was extracted, using a cir- 
cular region about 2' radius around the source centroid. 
The background spectrum was extracted from a sur- 
rounding annulus of radii 4 '.3 and 7 '.5, after removing 
a 2' circular area around a weak serendipitous source at 
a 20 oo = 12 h 39 m 39 s .2, 5 200 o = -5°26'13". Total exposure 
time is 1261 s. Publicly available response matrices, ap- 
propriate for the epoch of the observation, were employed 
and spectral fit were performed in the energy range 0.1- 
2 keV, where these matrices are best calibrated. A sim- 
ple power-law model with photoelectric absorption is an 
adequate description of the spectrum (x 2 = 9/17 dof), 
with best-fit parameters: N H = (2.1 ± 0.5) x 10 20 cm' 2 , 
r = 2.13±g;$, 0.1-2 keV flux ~ 1.26x 10' 11 erg cm' 2 s -1 . 
The intrinsic spectral index derived by ROSAT is indeed 
steeper than measured by BeppoSAX or ASCA, but by 
an amount much smaller than reported by Walter & Fink 
(1993). The soft excess above a T = 1.86 power-law is only 
20% in flux against the 270% reported by these authors. 
The different energy band-passes in which the spectra are 
taken, or slight residuals misalignment in the cross cal- 
ibration between ROSAT detectors and other missions 
(cf. Yaqoob et al. 1994; Iwasawa et al. 1998; Iwasawa et 
al. 1999) may easily account for this effect. Interestingly 
enough, no absorption edge is required by the ROSAT 
data despite the better effective area of the PSPC in com- 
parison to the LECS (the energy resolution is comparable 
at 0.5 keV). The upper limit on the optical depth of a 
0.77 keV absorption edge is 0.26. This might be sugges- 
tive of a long-term positive correlation between flux and 
warm absorber features, since the 0.1-2 keV flux is 60% in 
the ROSAT observation than in the BeppoSAX one. This 
aspect may be worth further investigation by future X-ray 
monitoring programs. 

6. Discussion 

The broadband (0.1-200 keV) X-ray spectrum of the 
Seyfert 1 galaxy NGC 4593 shows most of the typical 
features of its class. The continuum is well described by 
the superposition of a power-law primary component with 
r ~ 1.9 and a Compton reflection component. The spec- 
tral index is broadly consistent both with Ginga (Nandra 
k Pounds 1994) and ASCA (Nandra et al. 1997) mea- 
surements, and close to the typical values observed in 
Seyfert Is as a class (Nandra & Pounds 1994, Nandra et al. 
1997). The best- fit nominal values on R is consistent with 
a slab geometry for the reflecting matter. The reflection 
component is therefore likely to originate in an X-ray illu- 
minated relativistic accretion disk (AD), which subtends 
a ~ 27r solid angle from the nuclear source. This picture 
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is consistent with the properties of the iron emission line. 
Its centroid energy is well consistent with Ka fluorescence 
from neutral or mildly ionized iron (E c = 6.42±q }|). The 
line is moderately broad (a ~ 0.3). Fluorescence iron lines 
are expected to be produced, along with the Compton 
reflection continuum, in AD and broadened by the com- 
bination of gravitational and Doppler effects if the pho- 
tons undergo the effects of the gravitational potential of 
a supermassive black hole. The ~ 200 eV EW is around 
the expectation values if the AD matter has nearly solar 
abundances, given the measured value of R (Matt et al. 
1992). NGC 4593 does not suffer the iron overabundance 
problem that affects the relativistic iron line measure in 
some bright objects (among which the best studied case so 
far: MCG-6-30-15, Tanaka et al. 1995; see also Nandra et 
al. 1997) or to require a contribution from a narrow line 
component, which may originate in the molecular torus 
surrounding the nuclear environment in the unification 
scenario (Ghisellini et al. 1994; Krolik et al. 1994), and 
observed in several cases (Guainazzi et al. 1996; Weaver 
et al. 1997; George et al. 1998; Guainazzi et al. 1998a). 
The line properties can in principle constrain the location 
of the emitting region and the geometry of the AD. Fitting 
the iron line profile with a relativistic model allows us to 
constrain the inclination of the system (9 = 32° ±^), un- 
der the assumptions that the bulk of the iron line emitting 
region is neutral. 

The combination of intermediate X-ray continuum 
and iron emission line properties is hence perfectly con- 
sistent with the standard picture for the production of 
high-energy radiation in the nuclear region of radio-quiet 
nearby AGN. The primary continuum, with the "canon- 
ical" r = 1.9 power-law index, is isotropically produced 
in the neighborhood of the nuclear supermassive black 
hole and reprocessed via Compton down scattering and 
fluorescence by an optically thick and geometrically thin 
Shakura-Sunyacv disk (Shakura & Sunyaev 1973), whose 
extension is much larger than the typical size of the pri- 
mary continuum production region. The reprocessing mat- 
ter does not exhibit any substantial ionization or deviation 
from cosmic abundances. 

EXOSAT (Ghosh & Soundararajaperumal 1993; 
Santos-Lleo 1994) observed a strong and remarkably vari- 
able soft excess. Its intensity was between and 45% of 
the extrapolated high-energy flux in seven EXOSAT ob- 
servations, when soft and intermediate X-rays were mea- 
sured simultaneously. In EXOSAT data, it seems to exist 
no obvious correlation between the soft excess and the 
intermediate X-ray flux or spectral index (Santos-Lleo et 
al. 1994). Santos-Lleo et al. (1994) discuss the soft excess 
in terms of thermal emission from an accretion disc and 
conclude that the maximum temperature T max should be 
~ 210 cV. Such a disk blackbody component would be 
clearly detectable by the LECS (if not by the ASCA/SIS). 
If we add to the BeppoSAX "baseline" model the emission 
from a disc blackbody with temperature at the innermost 



disk radius T = 100 cV = 0.488T max (Pringle 1981), its 
luminosity is < 3 x 10 42 erg s _1 . A possible physical ex- 
planation for the disappearance of the soft excess may be a 
cooling of the disc temperature profile, which shifts T max 
out of the BeppoSAX sensitive bandpass. Alternatively, 
the lack of soft excess detection could be simply due to 
the improved continuum determination, made possible by 
the broadband BeppoSAX spectral coverage. The 2-6 keV 
spectrum, against which Santos Lleo et al. (1994) mea- 
sured the soft excess in EXOSAT data, was flatter than 
in our BeppoSAX baseline model (< Texosat — 1-73 >). 
If one fits the LECS/MECS data between 2 and 6 keV, 
one gets indeed T ~ 1.76. Extrapolating this power-law in 
the 0.1-2 keV BcppoSAX/LECS leaves a complex residual 
spectrum, with positive and negative wiggles, dominated 
by the Ovil absorption feature. The interpretation of this 
residual spectrum in the broad EXOSAT/LE filters is not 
totally unambiguous. We conclude that the claimed soft 
excess might be simply an artifact, due to extrapolating 
the 2-6 keV spectrum into the low energy filter EXOSAT 
energy bandpass and/or to fitting a complex ionized ab- 
sorber with the Galactic cold photoabsorption only. On 
the other hand, the lack of soft X-ray excess is in line 
with the lack of blue bump in this source (Santos-Lleo et 
al. 1994). 

Another new result emerging from the BeppoSAX ob- 
servation of NGC 4593 is the lack of any significant cutoff 
on the primary power-law, a lower limit on E cuto ff being 
~ 150 keV. This outcome is consistent with the measures 
of cutoff energies on individual sources available so far: 
NGC 4151 (70-270 keV, Zdziarski et al. 1995; Piro et 
al. 1998), IC4329A (240-900 keV Madejski et al. 1995), 
MCG-6-30-15 (100-390 keV; Guainazzi et al. 1999), and 
NGC 5548 (110-200 keV; Nicastro et al. 1999). 
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